Analytical approach to semiconductor Bloch equations 
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PACS 71.35.-y - Excitons and related phenomena 

Abstract. - Although semiconductor Bloch equations have been widely used for decades to 
address ultrafast optical phenomena in semiconductors, they have a few important drawbacks: (i) 
Coulomb terms between free electron-hole pairs require Hartree-Fock treatment which, in its usual 
form, preserves excitonic poles but loses biexcitonic resonances, (ii) Solving the resulting coupled 
differential equations imposes heavy numerics which completely hide the physics. This can be 
completely avoided if, instead of free electron-hole pairs, we use correlated pairs, i.e., excitons. 
Their interactions are easy to handle through the recently constructed composite-exciton many- 
body theory, which allows us to analytically obtain the time evolution of the polarization induced 
by a laser pulse. This polarization comes from Coulomb interactions between virtual excitons, 
but also from Coulomb-free fermion exchanges, which are dominant at large detuning. 



Introduction. — Interaction of semiconductors with 
ultrafast laser pulses has been widely studied for decades. 
Light-matter interaction leads to a large variety of non- 
linear effects of great interest in semiconductor technol- 
ogy. Probably the most fascinating ones are those in- 
duced by unabsorbed photons: while tliese photons could 
be thought of being completely inactive, they do act on 
semiconductors in a quite subtle manner through their 
coupling to virtual excitons. These virtual excitons inter- 
act with carriers, either electrons or holes already present 
in the sample, not only through the Coulomb potential, 
but mostly through the Pauli exclusion principle. The 
many-body theory for composite excitons we have recently 
constructed [1] makes this Pauli exclusion principle ap- 
pear by means of "Pauli scatterings" between excitons. 
They correspond to fermion exchange in the absence of 
fermion interaction, as nicely visualized by the so-called 
Shiva diagrams [2] . These Pauli scatterings in fact control 
all nonlinear optical effects in semiconductors. Among 
those induced by unabsorbed photons, we can cite the 
exciton optical Stark shift [3] , Faraday rotation [4] in pho- 
toexcited semiconductors, precession [5] and teleportation 
[6-8] of trapped electron spin, quantum computing [8-10] 
and classical optical computing [11] using trapped spins. 
Effects induced by these virtual excitons decrease with 
photon detuning but increase with photon number, i.e., 
pulse intensity. 



From a theoretical point of view, these effects face a ma- 
jor problem. As excitons play a key role, we must include 
electron-hole Coulomb interaction exactly, in order to pos- 
sibly reach their associated poles. However, many-body 
effects are known to be approached through perturbative 
expansion only. This is fully reasonable for interaction be- 
tween the exciton fermionic components, since excitons are 
neutral objects, so that their Coulomb interaction should 
induce small effects. Nevertheless, an exact treatment of 
the electron-hole interaction is required to get the exci- 
tonic poles. The clean separation between electron- hole 
ladder processes binding the exciton and those entering 
exciton-exciton interaction has stayed an open problem 
for decades [12]. One of the advantages of the composite- 
exciton many-body theory [1] is to provide such a separa- 
tion. 

In addition to the difficulty linked to the electron-hole 
interaction, excitons also have an intrinsic difficulty linked 
to their compositeness, which leads to a serious ambiguity 
when trying to pair up electrons and holes to form ex- 
citons. We have shown [1] that, as a direct consequence 
of the Pauli exclusion principle, it is impossible to prop- 
erly handle the interaction between excitons through an 
effective exciton Hamiltonian reading as Hx + Vxx ■ The 
potential Vxx, which describes interactions between ex- 
citons considered as elementary bosons, contains exciton- 
exciton scatterings. These are derived from the Coulomb 
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scatterings between elementary carriers through bosoniza- 
tion procedures which include some fermion exchanges. 
With the Hamiltonian not written as Hx + Vxxi it is 
not possible to use standard many-body techniques since 
they all rely on perturbative expansion in interaction po- 
tential. A totally different path had to be constructed to 
handle many-body effects with cxcitons. The many-body 
theory we proposed [1] is quite original as it uses an oper- 
ator algebra, instead of the scalar algebra used in Green's 
functions. The idea is to include the major difficulty of 
the problem, namely, the exact summation of all electron- 
hole processes giving rise to excitonic poles, in the very 
first line of the theory, by working with exciton opera- 
tors. These operators have an important drawback since 
excitons are composite bosons with commutation relations 
not as strict as for elementary particles. We can however 
deal with this difficulty by introducing "Pauli scatterings" 

A ■^^ which describe fermion exchanges in the absence 
of fermion interaction. They appear in 
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Since exciton operators b\ contain the tricky electron- 
hole part of the Coulomb interaction leading to exci- 
tonic poles in an exact way, the rest of Coulomb interac- 
tion can be treated perturbatively through exciton-exciton 

Coulomb scatterings ^ ■'^ . They appear in 



[h,b\] = e,b\ + v^ 



t nt 



(3) 
(4) 



differential equations. This is definitely a great advan- 
tage because the physical meaning of the various contri- 
butions to the polarization then appear in a transparent 
way. We find, once more, that one part of the polarization 
comes from Coulomb interaction between excitons while 
the other part comes from the Pauli exclusion principle be- 
tween the exciton fcrmionic components, which, as usual, 
is dominant for large photon detuning. 

Proposed approach. — We calculate the polariza- 
tion induced by a single laser pulse starting from the non- 
excited semiconductor. Theoretical extensions of this ap- 
proach to semiconductors already containing carriers or 
to experiments with more than one type of photons, such 
as in four- wave mixing, can be done along the same line. 
These will be presented elsewhere, as the goal of this let- 
ter is to present the main steps of the new procedure we 
propose. 

The time evolution of the exciton polarization can be 
written in terms of (V't|-Bp|V't) = (^p)ti where \ipt) is the 
semiconductor state resulting from its coupling Wt to the 
photon field. This mean value obeys 



{[H,, + Wt,Bl])t 



(7) 



is the semiconductor Hamiltonian in the absence of 
photon, while Wt = Ut + U}, with Ut = E» ^iWB^, the 
time dependence of the Rabi coupling fli{t) to exciton i 
following the laser pulse. Since for composite-boson oper- 
ators [B^,bI] — 0, so that [Uj ,Bp\ = 0, the commutation 
relations in Eqs.(l,3) readily turn Eq.(7) into 

^§l-E,^ {Bl)t^np{t) 

-Y.a,{tm,)t + {v^)t. (8) 



These Pauli and Coulomb scatterings allow us to calculate 
any physical quantity dealing with excitons in an easy way. 

Before using this new many-body theory to calculate 
the time evolution of the polarization induced by a laser 
pulse, we wish to emphasize that the semiconductor Bloch 
equations (SBE) [13-20] use free pairs ajj^^k^ instead of 

excitons Bj. This just corresponds to a basis change, since 
these operators are linked by 



Bl= E «L&L(k/.,keN) 



(5) 



(6) 



Difficulties with SBE actually come from using these 
free-pair operators in problems in which excitons play a 
key role. We here show, by working with exciton oper- 
ators in the form {B^)t instead of (Ok&Lk)*' 
calculate such an expectation value analytically, thereby 
avoiding the problem of solving numerically the coupled 



Beside a bare linear term flp{t) in photon coupling, the 
time evolution of {B^)t contains contributions coming 
from the composite nature of exciton p through Dip, and 
from Coulomb interaction with exciton p through . 

If we only keep resonant terms, for a laser field with 
frequency lo introduced adiabatically from t = —oo with 
a 1/e risetimc, the time dependence of the Rabi coupling 
ili{t) reads rj^e'^'e'*, so that, for a non-excited semicon- 
ductor, \^t=-oci) = \v), the linear term rip(t) in Eq.(8) 
gives the free part of (i?^)t as ripe*'^*e'^*/ {uj — Ep — ie). 

However, even if adiabatic establishment is considered 
in most papers, it does not seem to us the proper way to 
describe ultrashort laser pulses with steep rise time and 
duration time Tp long compared to photon period 
We find it more appropriate to take ^i{t) as n^e*"* during 
the pulse duration, < t < Tp and zero otherwise, with 
iV't) = K') for t < 0+ (sudden approximation, valid for 
steep rise time). The free part of (i3p)t, which now cancels 
for t = instead oi t = —oo, reads as 



HBlYr^e- 



(i; 
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where ujp ^ lu — Ep is the photon dctunmg with respect to 
excitonp, while A^^^(a;) = (e''^* - l)/x, so that A^^'(O) = 

it. Note that, as A[^\x) also reads 2i'Ke^^*^'^St{x), where 
St{x) = {ttx)~^ sin{xt/2) is a "delta function" of width 
2/t, this free part barely leads to the Fermi golden rule. 

In addition to its physical relevance, a sudden rise 
for iliit) allows us to eliminate t from the Hamiltonian 
through the so-called "rotating frame" unitary transfor- 
mation that we can take [21] as Zt = exp{—iujt^a^ai^). 
This gives Zf^^a^^Zf = e*'^*aj^, while stays unchanged. 
The effective Hamiltonian H = Z^^HtZf — iZ^^Zt, which 
rules the time evolution of j-i/jj) ~ Z^^ltpt), reduces to 
Hsc + W, where Hgc is just Hgc with all electron energies 
shifted by —to, which amounts to replace Ei by Ei —uj, with 
the coupling W reading asU + W where U — ^iBi is 
now t-indepcndent. The time evolution of 



{MBl\^t) = {{B%^e-^-\Bl), 



(10) 



follows from Eq.(8), with —Ep replaced by ujp and ^i{t) 
by 

The linear term Q.p leads to {{B^p)Yt'"^, given in Eq.(9). 
To get Pauli and Coulomb contributions coming from 
{{Dip))t and {{V^))t, we note that, since H is now time- 



independent, we then have \ipt) 



-iHt 



v), in which we 



can use the integral representation of the exponential, 



iHt 



dx 



-2iTr) 



^ — i{x-^iO-i^)t 

+ iO+-H 



(11) 



and decouple interaction with photons from interactions 
between carriers, through 



1 



1 



a — H a — Hsc a — H a — Hsc 

= y^f^^wX (12) 

„=o ^ " ^ ^ "■~ ^''c 

As W changes the number of pairs by one while Dip 
conserves this number. ((Dip))t only has even terms in 
photon coupling, while {{Vj^))t only has odd terms since 
Vj^ creates a pair. Therefore, as expected, the polarization 
given in Eq.(8) has odd terms only in Rabi coupling Qi. 

Eq.(12) generates the polarization as an expansion in 
photon interaction W. Such an expansion is a priori valid 
for small laser intensity. We wish to stress that this limi- 
tation is fully consistent with only considering {{B))t but 
not ((_B"))f with n > 2. In the case of large photon field 
for which such a W expansion is not valid, the relevant 
operators are no more exciton, but polariton operators. 
In a near future, we will address large photon fields us- 
ing our work on interacting polaritons [22,23], along ideas 
similar to the ones we propose here. However, polaritons 
being far more complex composite bosons than excitons, 
this extension is at the present time beyond our scope. 



Pauli part of {{Bj^))t. — As shown in the appendix, 
the Pauli part of {{B^))t, which comes from the second 
term of Eq.(8), has a third order contribution in photon 
coupling which reads 

((i?t))f-" ^ 2e-'^-'J2^^^^^k>^^jkpM'\^P.-^j,^k) , 

ijk 

(13) 



where 2Xijkp 
tions are linked by 



A (^ p')+(*-j) 



while the aI"'' func 



A(n+1)/ \ 
t '{Xi,X2,---,Xn+l) = 

A^"^(a;i + .T„+i,a::2,---,a;„) - A^"^(a;i, a;2, • • • , a;„) . (14) 

Coulomb part of {{Bjj))t. — The Coulomb part of 
{{Bj,))t, comes from the last term of Eq.(8). In contrast 
with {{Dip))t, which only reads in terms of exciton energies 
and wave functions, the exact calculation of {{Vj))t at 
third order in photon coupling requires the knowledge of 
the whole two-pair eigenstate spectrum. We can however 
approximate this Coulomb part in two limits: 

(i) For materials having a well-separated biexciton, the 
Coulomb part of the polarization for photons close to the 
biexciton resonance is found to read 



{{Bl))Y"' 



ijk 



{XX) 
ijkp 



X A^^Up^Exx -2lj + Ljj,-u)j,ujk) ; 



(15) 



f (^"^' can be seen as a Coulomb interaction between ex- 
citons mediated by the molecular biexciton resonance. 



Y^{v\B,B,\XX){XX\BlBi\v)^{^l) , (16) 



where |^^) and Exx are the molecular biexciton state 
and energy, these being obtained through numerics only. 

(ii) We can also expand the Coulomb part of the po- 
larization in Coulomb scattering divided by detuning, 
through 



1 



Hsc 



B 



1 



1 



a — J^sc " \ " a — Hsc J a ^ Hsc + 
As shown in the appendix, we then find 

((Bt))Coulomb 



(17) 



ijk 



X A[^\cup, 



i, ~^],^k)) , 



(18) 



where ^ is the mixed direct-exchange Coulomb scattering 
standard for time evolution [1] , 

24-fcp = [cOp)-rOp)] + (*-j), (19) 

the "in" Coulomb exchange part of this scattering being 

equal to E™,„a(^ m) ^ (mp)- 
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Discussion. — The above derivation shows that 
{{Bjj))t contains two cubic terms in photon coupling given 
in Eq.(13) and Eq.(18) — or Eq.(15) when the molecular 
biexciton plays a key role. A term comes from fermion ex- 
change with exciton p without Coulomb process, the other 
from Coulomb interactions with exciton p with or without 
fermion exchange. They basically have the same structure 
except that, Coulomb scatterings ^ijkp being energylike 
quantities while Pauli scatterings Xijkp are dimensionless, 
the Coulomb part of the polarization has one more dctun- 
ing in the denominator than the Pauli part: indeed, AJ 

has LUjUJk, while Aj^'' has ujiUJjOJk- Consequently, as usual 
for nonlinear optical effects, the large detuning behavior of 
{{Bjj))t is fully controlled by the Pauli exclusion principle 
with exciton p in the absence of Coulomb interaction. 

Furthermore, we note that the coupling Qi to exciton i 
depends on i through (r = 0\i/i), where (r|i/i) is the i ex- 
citon relative motion wave function. This makes photons 
dominantly coupled to s excitons. In addition, ((i?j,))t de- 
creases when photon detuning increases. Therefore, size- 
able {{Bj^))t comes from photon detuning to the Is exciton 
ground state (labelled 0) not too large. For such photons, 
the uJiOjjUJk denominator of AJ in Eq.(18) leads us to only 
keep {i,j,k) = in the sum. Similarly, the ujjUJi; denomi- 
nator of Ap-* in Eq.(13) leads us to only keep {j,k) = 0. 
The cubic term of {Bj^)t then appears as 



(4) 



t\(3) 



(3) 



(20) 



where r^oAooop = X^i ^lAioop, the summation over i being 
performed through closure relation. The above result is 
valid for uiq = uj — Eq small compared to uii^o but still large 
enough to possibly keep one term only in the ((-8^,))^°^'°™'^ 

expansion in ^/wq. For p being a bound state, Aooop — 
Xoiax/L)^ and |ooop = CoRxiax/L)^, where ax and Rx 
are the exciton Bohr radius and Rydberg energy, D the 
space dimension, L the sample size and Xd, numerical 
prefactors of the order of 1. Analytical values of these 
prefactors for p = can be found in refs. [1, 22]. 

For t small, i.e., for ujt ^ 1, we do have A^"-* ~ (it)"-/n, 
so that (Bj,)^'' rises as 



(3) 



2A 



OOOp- 



{^tf 



Cooop 



(21) 



More generally, the recursion relation between the 
a|"'''s gives the t dependence of the Pauli part of (i?^)t as 

u^^ [a|i^ (c^p + c^o) + a|^) [ujp - uj^) - 2A W {ujp)\ , (22) 
while for the Coulomb part, we find 

^t'^' {^Pi -^0, -Wo, Wo) = Wq"^ 3A^^^ {up - Wo) 



-f A^^^ {ujp + ujq) ~ a|^^ {up - 2a;o) - SA^^^ (wp)] . (23) 

where Aj^' (w) is essentially a delta function in the large t 
limit. Fig.l shows the real parts of a|"'' for n = (1,3,4), 
when p — Q. 

State of the art. — Polarization is usually derived in 
the free-pair basis by calculating = (ak&-k)t, through 
-«^Pk = {[{Hsc + VFt), Ok^-k])*- In this commutator, 
the Coulomb part of Hsc generates expectation values like 
{a^a^aa)t, {b^b^bb)t and {a^b^ba)t- These are then cut by 
Hartree-Fock (HF) procedure, as necessary to get a closed 
set of differential equations between pk and carrier densi- 
ties, fl^^ = (aj^.ak)t and j^'' = (&k&k)t- These equations 
are then solved numerically using a 4*'' order Runge-Kutta 
scheme. 

This SEE procedure is presented in details in ref. [13]. 
It however has some unpleasant drawbacks, such as di- 
vergence in carrier-carrier scattering rates coming from 
HF approximation, which has to be cut in a "controlled" 
way. Standard HF approximation allows one to recover ex- 
citonic poles [13], but higher-order Coulomb correlations 
are out of reach. These, including interband electron-hole 
exchange, can be taken into account, in the case of equi- 
librium, through the Bethe-Salpeter equation [24-28], or 
through nonequilibrium Green functions [29,30]. Both 
techniques are numerically quite demanding [28,31]. It 
has recently been found that single-exciton SBE based on 
time-dependent density functional theory (TDDFT) pro- 
vide a much simpler approach [32]. However, they can- 
not yet account for two-exciton correlations in the single- 
exciton polarization. In particular, dark excitons with 
spin (±2) are not reached [16] while, evidently, such dark 
excitons are produced by carrier exchange between two 
bright excitons [33,34]. Furthermore, the biexciton bind- 
ing energy cannot be reached from single-exciton SBE. 
Many important dynamical effects are thus completely 
lost in these equations. Possible improvement can come 
from dynamics-controlled truncation (DCT) at the two- 
exciton level in order to obtain biexcitonic correlation 
effects [13,17,18]. This DCT method is similar to the 
BBGKY (Bogoliubov-Born-Green-Kirkwood-Yvon) hier- 
archy in statistical mechanics for the time evolution of a 
many-particle system. In any case, since the obtained re- 
sults rely on heavy numerics, their physical understanding 
is difficult to grasp because these equations are supposed 
to contain "everything" at once, so that all effects are 
mixed in the final results, given through numerical curves. 

The procedure presented here is far more transparent 
since it allows us to find analytical expressions for the t 
dependence of the excitonic polarization, up to any arbi- 
trary order in Coulomb processes and photon field — this 
photon field expansion being fully consistent with con- 
sidering one-pair operators instead of polariton operators. 
Our new procedure has the huge advantage to trace back 
all dynamical effects to these expansion terms. It is in 
particular possible to distinguish effects arising from Pauli 
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Fig. 1: (a): Real part of t''^ A[^\u}o) , defined in Eq.(9), as a function of uot. (b) and (c): Real parts of t ^A[^\ljo, —loo,ujo) and 
t~^A['^\u!o, —too, —iuojUJo) defined in Eqs.(22,23), A^"' having the same parity as n. These curves are more and more picked 
around = when t or n increases. 



blocking, from effects arising from Coulomb interaction. 

Conclusion. — We here develop an alternative ap- 
proach to semiconductor Bloch equations which allows us 
to calculate the time evolution of the polarization induced 
by a laser pulse analytically. This approach completely 
avoids the heavy numerics associated with solving these 
equations as well as the spurious singularities originating 
from their necessary truncations. 

Such an analytical solution has been made possible be- 
cause we here use correlated pairs, i.e., excitons, instead of 
free pairs, their interaction, responsible for the polariza- 
tion, being handled through the new composite-exciton 
many-body theory [1]. This analytical approach makes 
the two physical channels producing the polarization com- 
pletely transparent: fermion exchange and fermion inter- 
action, the former being dominant at large photon detun- 
ing. Our approach also deals with the biexcitonic reso- 
nance in an easy way, this resonance being uneasy to reach 
within the usual semiconductor Bloch equation formalism. 
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Appendix. — For readers interested in technical as- 
pects of this approach, we now give the main steps leading 
to Eqs.(13,15,18). 

(a) The Pauli part of {{Bj^))t comes from {ipt\Dip\ipt) . 
By using Eqs.(10,ll), this mean value can be written as 



2i7r 



-f-oc 



dx 



J{x'-iO+)t -i{x+iO+)t 



(24) 



where, since Dip\v) = while Hsc\v) = 0, we do have 



^ +00 

Vip{a' , a) ^ -J— ^ {4'n'{a')\Dip\(j)n{a)) 



(25) 



l0Ti(a)) 



W 



(26) 



^ a — Hsc / d ~ 

Due to Eqs.(l,2), {v\BjD,pBl\ v) is simply equal to 2Xijkp- 
This gives the quadratic term in photon coupling of 
Dip{a' , a), which comes from n = n' = Q in Eq.(25), as 



a' a ^-^ (q! + + ^fc) 

jk 



A. 



ijkp 



(27) 



(2) 

{{Dip))\ follows from the above equations inserted in 
Eq.(24). Eq.(8) then leads to Eq.(13) for the cubic con- 
tribution to {{Bj,))t induced by Pauli blocking. Higher 
order contributions in photon coupling, which come from 
{n,n') ^ (0,0) terms in Eq.(25), can be obtained in the 
same way. 

(b) The Coulomb part of {{Bj^))t comes from {{Vj))t 
which reads as {{Dip))t in Eqs.(24,25), Dip being simply 
replaced by . Since creates a pair while Vj\v) = 0, 
the lowest order term in photon coupling, obtained for 
n = and n' = 1 instead of n = n' = as for Eq.(28), 
then gives 



Vf (a',a) = ^^- 

^ n' n ^— ' I n' 



a' a ^—^ (q! + + 

ijk 



Cypfe(a') . (28) 



Eq.(4) allows us to write 



Cijpk{a') = {v\BiBj — 



a' - 



'V^Bl\v) 



m.n 



(29) 



The exact calculation of 



My„„(a') = {v\B,Bj — 

rt' 



1 



BlBi\v) , (30) 



a' ~H,_ 

requires the knowledge of the whole two-pair eigenstate 
spectrum. 

(i) For materials with well-separated molecular biexci- 
ton, this particular two-pair eigenstate controls Mijmn{a') 
for photons close to the biexciton resonance. When in- 
serted into Eq.(30), the \XX) biexciton state gives 



ijnin 



(«') 



{v\B,B,\XX){XX\BlBi\ 
a' + 2lj- Exx 



(31) 
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Eq.(29) then leads to 



Cijpk (fl ) 



■ 2lu — Exx 



(32) 



with ^ijfe^'* given by Eq.(16). Inserting the above result 
into Eq.(28), we get ((V^))t through an equation similar 
to Eq.(24), from which it is easy to find the Coulomb part 
of the polarization given in Eq.(15). 

(ii) Mijmnia') can also be calculated by iteration. In- 
deed, Eq.(17) leads to 



(a' + UJ-i + LOj)Mijjnn{a') = 



(«') 



(33) 



the first term in tlie above equation being equal to 
- -^(im) + (™ ^ due to Eqs.(l,2). By 
keeping this first term only, we get the first term of the 

tuning. Eq.(29) then gives 



-iijmn(a') expansion in Coulomb scattering divided by de- 



(34) 



with (,ijkp defined in Eq.(19). The third-order contribution 
in photon coupling coming from Coulomb interaction with 
exciton p, as given in Eq.(18), then follows from Eq.(28). 
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